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O ■ ABSTRACT 



After a short introduction on the predicted cosmic neutrino background in the 
universe, we review some of the cosmological bounds related to neutrinos. In 
particular we show how the recent data on the anisotropics of the cosmic mi- 
■ crowave background radiation from the WMAP satellite, combined with other 

(SJ , experimental results, can constrain the radiation content of the universal energy 

density and the sum of neutrino masses. 

Cf) , 1. Introduction 

m 

Neutrinos are the second most abundant particles in the universe, with a number 
density comparable to that of relic photons. Therefore, the relic neutrino sea influ- 
ences various cosmological situations, playing an important role that has been dis- 
cussed in many papers, both to solve some cosmological problems and to put bounds 
on non-standard neutrino properties. These cosmological limits are alternative to the 
o ■ data from terrestrial neutrino experiments, in some cases the only available. 

Here we review recent bounds from cosmology on the number of neutrinos (content 
of relativistic particles) and on the neutrino mass, in particular those that appeared 
after the release of the first year data of the Wilkinson Microwave Anisotropy Probe 
(WMAP) EJ. There exist many other topics not discussed here, such as the role of 
massive neutrinos in leptogenesis, that can be found in the extensive review ^ and 
in the more recent ones 1"™ . 

2. Standard relic neutrinos 

In the early universe, at temperatures above MeV, the neutrinos were in thermal 
equilibrium through the standard weak interactions with other particles. Thus the 
distribution of neutrino momenta was a Fermi-Dirac one, 



UJp) 



exp ( ) + 1 



where the energy E ~ p for neutrino masses much smaller than MeV. The the 
chemical potentials of the different neutrino flavours, which could have been non-zero 

a Invited talk at the X Int. Workshop on Neutrino Telescopes, Venice, March 11-14, 2003. 



if an asymmetry between the number of neutrinos and antineutrinos was previously 
created. We neglect this possibility in this section, but it will be considered later. 

As the universe cools, the weak interaction rate T v falls below the expansion 
rate H and the neutrinos decouple from the rest of the plasma. A rough but quite 
good estimate can be found equating the average value of T u = (cr u n u ) and H = 



8np/3Mp, which gives T dec ~ O(MeV). A better calculation gives T dec ~ 2 — 3 
MeV, where the higher (lower) decoupling temperature corresponds to v^ T (z/ e ). After 
decoupling, the collisionless neutrinos expand freely, keeping their phase-space density 
as in Eq. JIJ), since both momentum and temperature redshift identically. Finally, at 
temperatures of order the electron mass, the electron-positron pairs annihilate and 
transfer their entropy into photons but not into the decoupled neutrinos, causing the 
well-known difference between the temperatures of relic photons and relic neutrinos, 
T 7 /T,= (ll/4) 1/3 (seee.g.EJ). 

Note that the relic neutrino phase-space density in Eq. (JTJ) is that of relativistic 
particles (or radiation) in equilibrium, valid also for massive neutrinos when m v ^> T v . 
It is easy to calculate the number and energy densities of relic neutrinos at any 
temperature T = T 7 below MeV. The number density per flavour is fixed by the 
value of the temperature, 

n-^n=^T* (2) 

while the energy density for massive neutrinos should be in principle calculated nu- 
merically, with two well-defined analytical limits, 

pM» « T v ) = — [-) T 

Pu(m u > T v ) = m u n u (3) 

Thus we see that the contribution of massive neutrinos to the energy density in the 
non-relativistic limit is a function of the mass (or the sum of masses if all neutrino 
states have mj ^> T v ). 

Relic neutrinos contribute to the relativistic energy density of the universe while 
m u T u , fixing the expansion rate during the production of primordial abundances 
of light elements (Big Bang Nucleosynthesis or BBN). BBN is the last period of the 
universe sensitive to neutrino flavour, since electron neutrinos and antineutrinos play 
a direct role in the freezing of the neutron-to-proton ratio. Later cosmological stages, 
such as the formation of the angular fluctuations of the cosmic microwave background 
(CMB) and the formation of large-scale structures (LSS), are also sensitive to rela- 
tivistic relic neutrinos but blind to neutrino flavour. In particular the small-scale 
power of the cosmic density fluctuations is modified by neutrino masses of order eV 



and smaller. 



3. Number of relativistic neutrino species 

The energy density of the universe stored in relativistic species, p r , is customarily 
given in terms of the so-called effective number of relativistic neutrino species N e g 
(see^ for a review and references), through the relation 



Pr = P 7 + Pu + Px 
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1 + -(n) v 
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P7 , ( 4 ) 



where p 7 is the energy density of photons, whose value today is known from the 
measurement of the CMB temperature. Eq. (jH) can be also written as 



N* = ^ - ■ (5) 




where p% denotes the energy density of a single species of massless neutrino with 
an equilibrium Fermi-Dirac distribution with zero chemical potential, and p° is the 
photon energy density in the approximation of instantaneous neutrino decoupling. 
The normalization of iV e fT is such that it gives N e g = 3 in the standard case of three 
flavours of massless neutrinos, again in the limit of instantaneous decoupling. In 
principle N e g includes, in addition to the standard neutrinos, a potential contribution 
p x from other relativistic relics such as majorons or sterile neutrinos. 

It turns out that even in the standard case of three neutrino flavours the effective 
number of relativistic neutrino species is not exactly 3. The decoupling of neutrinos 
from the rest of the primordial plasma occurs not far from the e^ annihilations, and 
accurate calculations^ have shown that neutrinos are still slightly interacting, thus 
sharing a small part of the entropy release. This causes a momentum dependent 
distortion in the neutrino spectra from the equilibrium Fermi-Dirac behavior and a 
slightly smaller T/T u ratio. Both effects lead to a value of N e g = 3.034. A further, 
though smaller, effect on T/T v is induced by finite temperature Quantum Electrody- 
namics (QED) corrections to the electromagnetic plasma. A recent combined study 
of the incomplete neutrino decoupling and QED corrections concluded that the total 
effect corresponds to N e g = 3.0395 ~ 3.04 . Therefore we define the extra energy 
density in radiation form as 

AN cS = N cS - 3.04 . (6) 

The standard value of N e g corresponds to the case of massless or very light neutrinos, 
i.e. those with masses much smaller than 1 eV. More massive neutrinos affect the 
late evolution of the universe in a way that can not be parametrized with a AN e g, as 
described in the next section. 
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Figure 1: The ler, 2a and 3<r contours in the 7710 — AiVeg plane from a combination of BBN and 
WMAP data. Here uj b = Q B h 2 and r? 10 = 10 10 ^ B ~ 274 w B . Figure from ref.ES. 

The value of AiVeff is constrained at the BBN epoch from the comparison of 
theoretical predictions and experimental data on the primordial abundances of light 
elements, which also depend on the baryon-to-photon ratio 77 = ns/n^ (or baryon 
density). The main effect of a non-zero AN e g is to modify the Hubble expansion rate 
through its contribution to the total energy density H = ^8np/3Mp. This in turn 
changes the freezing temperature of the neutron-to-proton ratio, therefore producing 
a differen t abundance of 4 He. Typically, the BBN bounds are of order AN e s < 0.4 — 1 
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The BBN bounds on AiV e g- have been recently reanalyzed using as input the value 
of the baryon density derived from WMAP data^^ ?7cmb — 6.14 ± 0.25. Hannestad 
^ fin ds the constraint N eS = 2.6±g;| (95% CL), in agreement with the results of refs. 
| 14 | 15 | r-Q^ allowed regions are shown in Fig. Q taken from ref. The bounds 
suggest that the preferred value of AiV e g is below zero, i.e. inconsistent with the 
standard case. However, it is probably premature to consider this a robust indication 
of non-standard physics, because the observational data on primordial 4 He is still 
dominated by systematics. It is also interesting that, since the observed primordial 
abundance of deuterium agrees well with the value predicted from t/qmb, one can 
now use D to obtain a limit on N e g. The authors of found, for an averaged D 
abundance, the limits 1.22 < N cS < 5.25 (95% CL), in agreement with the standard 
prediction. In the future, improved bounds can be obtained for more precise values 
of ?7cmb and the primordial abundance of D. 

It is important to note that the BBN bounds on AiV e ff change if the v e or v e 
spectrum show deviations from the standard Fermi-Dirac form with [i Ue = 0. For 
details, see for instance the recent discussion in^^. 
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Figure 2: The two-dimensional confidence limits on (AA^ff , h), based on CMB and LSS data, at the 
1-cT (dark blue) and 2-a (light blue) levels. The superimposed yellow stripe shows the HST result. 
Figure from ref. . 



Independent bounds on the radiation content of the universe at a later epoch 
can be extracted from the analysis of the power spectrum of CMB anisotropies. An 
enhanced contribution of relativistic particles delays the epoch of matter-radiation 
equality, which in turn increases the early integrated Sachs- Wolfe effect. Basically 
this leads to more power ar ound the sca le of the first CMB peak. Previous analyses 
found weak bounds on ^jv eff II2IIBII113 ) that can be significantly improved by adding 
priors on the age of the universe or by including supernovae and LSS data^. Before 
the WMAP data, from a combination of CMB and LSS data from the PSCz redshift 
survey, ref. found the allowed range N cS = 6^4 5 (95% CL). 

Th e inclu sion of the WMAP data has significantly tightened the bounds, as shown 
in refs. 1 1 Assuming a minimal set of cosmological parameters and a flat universe, 
the limit is now AN eS = 0.5±l;f (95% CL), a value found using data from WMAP and 
other CMB experiments, as well as LSS data derived from the 2dF Galaxy Redshift 
Survey (2dFGRS). In particular, a parameter degeneracy exists between AN e s and 
the present value of the Hubble constant h = i/o/(100 km s _1 Mpc -1 ), which is 
par tially removed including the results of the Hubble Space Telescope Key Project 
(HST), h = 0.72 ± 0.08 at la. This degeneracy is shown in Fig. El taken from 
ref . ^21 (see also^ 1 ). In the case of a non-flat universe 123^ the corresponding limit 
is modified to AA" eff = l.ltij. 

Although the recent CMB+LSS bounds on AN e s are not as restrictive as those 
from BBN, they will be significantly improved in the future with the next measure- 
ments of the CMB temperature and polarization spectra by W MAP (an d those of 
PLANCK in the near future), as shown in several forecast analyses 122121123 An anal- 



ysis for models without dark energy contribution estimated that it would be pos- 
sible to measure AN e g ~ 0.04, the standard value calculated with non-instantaneous 
neutrino decou plin g. However, this is probably too optimistic, and a more recent 
estimate in ref.123 gives AN e ff ~ 0.2 with the inclusion of PLANCK data. 

Many extensions of the Standard Model of particle physics predict additional 
relativistic degrees of freedom that will contribute to AiV e fj. There exist models with 
4 neutrinos which include an additional sterile neutrino in order to explain the third 
experimental indication of ne utrino oscillations (the LSND results). It was shown in 
many studies (see for instance l ^ l ^ l) that all four neutrino models, both of 2+2 and 
3+1 type, lead to a full thermalization of the sterile neutrino before BBN, and thus 
to AiV e fT ~ 1, a value disfavoured in the standard minimal model of BBN. Moreover, 
in these models there exists at least one neutrino state with mass of order 1 eV, and 
some of the regions of parameters are also disfavoured by the cosmological upper limit 
on m ui as shown in the next section. 

It is also possible that the relativistic degrees of freedom at the BBN and CMB 
epochs differ, for instance because of pa rticle decays which increase the photon tem- 
perature relative to the neutrino one'^'. In some situations AiVpfr can b e effec tively 
negative at BBN, such as the case of a distortion in the v e or v e spectra \ 01 i °^ t or a 

F7*7I 

very low reheating scenario 1 ^. 

A non-standard case that has been considered many times in the past is the 
existence of relic neutrino asymmetries, namely when the number of neutrinos and 
antineutrinos of the same flavour is significantly different. These so-called degenerate 
neutrinos are described by a dimensionless chemical potential £ a = fi Ua /T, and it has 
been shown that the neutrino energy density always increases for any value £ a ^ 



AiV cff = ^E 




(7) 



Interestingly, some combinations of pairs (£ e , £ M)T ) could still produce the primordial 
abundances of light elements for a larger baryon asymmetry, in the so-called degen- 
erate BBN scenario . At the same time the weaker CMB bounds on £„ are flavour 
blind MB. However, it was recently shown that for neutrino oscillation parameters 
in the regions favoured by atmospheric and solar neutrino data flavour e quilibrium 
between all active neutrino species is established well before the BBN epoch 1HHISZISH1 
Thus the stringent BBN bounds on £ e apply to all flavours, so that the contribution 
of a potential relic neutrino asymmetry to AA e g is limited to very low values. In 
addition, these bounds fix the value of the neutrino number density down to 1%, an 
important input to find an upper limit on m u . 



4. Massive relic neutrinos 



The role of neutrinos as the dark matter (DM) particles has been widely discussed 



since the early 1970s. There are two facts in favour of neutrinos as DM: they definitely 
exist and an eV neutrino mass produces a contribution of order unity to the present 
energy density of the universe. This can be easily showed calculating the neutrino 
fraction of the critical density, Q u = p v j p c . For values of neutrino masses much larger 
than the present cosmic temperature (T 7 ~ T v pa eV), one finds 

From this relation, one gets the traditional upper limit on neutrino mas ses, of the 
order some tens of eV, using the very conservative bound Cl u < 1 (see 1*™ for a 
list of references). As discussed below, much tighter constraints are found from LSS 
observations. 

The evolution of cosmological perturbation is sensitive to the energy density of 
massive neutrinos. The neutrino background erases the density contrasts on wave- 
lengths smaller than a mass-dependent free-streaming scale. Thus massive neutrinos 
produce a damping of the total density fluctuations on small scales, characteristic of 
the so-called hot dark matter (HDM) particles. In a universe dominated by HDM, 
the formation of structures follows a top-down scenario, where large objects such as 
superclusters of galaxies form first, while smaller structures like clusters and galaxies 
form via a fragmentation process. 

In the late 1970s and early 1980s, neutrinos looked like a good DM candidate. 
But by the mid-1980s it was realized that HDM models were not able to explain 
some of the observed structures, among others the fact that galaxies seem older than 
clusters. The attention then turned to cold dark matter (CDM) candidates, i.e. 
particles which were non-relativistic at the epoch when the universe became matter- 
dominated, which provide excellent agreement with observations. Still in the mid- 
1990s it appeared that a small mixture of HDM in a universe dominated by CDM 
fitted better the observational data than a pure CDM model. However, within the 
presently favoured ACDM model (dominated by dark energy) there is no need for a 
significant contribution of HDM b . 

Nowadays we know that neutrinos can not constitute the main DM component, 
and the main question now is to find how large the neutrino contribution (or HDM) 
can be. For values of the neutrino masses of order eV, the free-streaming effect 
can be detectable in the linear power spectrum of LSS, reconstructed from red- 
shift survey experiments. This effect is shown in Fig. |2J where the suppression 
of the mass power spectrum P(k) caused by massive neutrinos was calculated for 
two cases with the same J2 m i, an d compared with the analytical approximation^^ 
AP(k)/P(k) ~ — 8Q U /Q m , usually quoted in the literature. Thus cosmology offers a 
unique opportunity to measure the absolute value of the neutrino masses. 



3 For a historical review of neutrinos as DM, see for instance 



39 1 




1 1 1 1 ^ 1 1 

0.5 1 1.5 2 2.5 3 

Em, 

Figure 3: Ratio between the power spectrum at k = 0.5ft. Mpc -1 with and without massive neutrinos, 
as a function of the total mass £ m„ in eV when the number of massive neutrinos is either three or 
one. The other cosmological parameters are kept fixed and the universe is assumed to be fiat. The 
analytical approximation often quoted in the literature is also shown (blue dashed). 

The free-streaming effect is more important than the direct effect of the contribu- 
tion of massive neutrinos to the energy density of the universe: different values of the 
neutrino density fraction Q v have to be compensated by different values of the other 
components (cold dark matter, baryons, dark energy or eventually the curvature pa- 
rameter). These changes in the cosmic balance can modify some characteristic times 
and scales in the history of the universe, like the time of equality between matter 
and radiation, or the size of the Hubble radius at photon decoupling. Neutrino mass, 
however, influences only slightly the spectrum of CMB anisotropics. Nevertheless, it 
is crucial to combine CMB and LSS observations in order to measure the neutrino 
mass, because CMB data give independent constraints on the cosmological parame- 
ters, and partially removes the parameter degeneracy that would arise in an analysis 
of the matter power spectrum only. 

These effects of neutrino masses probed by CMB and LSS depend, in first ap- 
proximation, on the sum of the masses m v if all neutrino states have the same 
number density, that of Eq. (j2J). Before WMAP data, the best bound on J2 m u was 
obtained using the power spectrum of galaxy clustering from the 2dFGRS and adding 
other independent cosmological constraints through priors on the other parameters 
22. The upper limit found was £ m v < 2.2 eV (95% CL). 

At present, the most stringent upper limit on the sum of neutrino masses was 
reported by the WMAP collaboration using the first-year WMAP data, CMB 



data from other experiments (ACBAR and CBI), the LSS data from 2dFGRS and the 
matter power spectrum on small scales inferred from the Lyman a forest. The limit 
is J2 m u < 0.7 eV (95% CL), approximately a factor three better than the previous 
2dFGRS bound. 

Two other recent analyses EMU have calculated the upper bound on m u within 
a more conservative approach. In particular the Lyman a data was not used, since 
it has been recently pointed out that the errors were probably underestimated ^31 
Hannestad also used a completely free bias parameter (which relates the matter and 
the galaxy correlation spectra), emphasizing the importance of the HST prior on h 
in constraining the neutrino mass. The upper bound found is 

J2m„ < 1.0 eV (95% CL) (9) 

which should be considered conservative, but still dependent on the set of cosmological 
parameters (a minimal set) and the assumed priors (see the discussion on 1 1 

It should be noted that this cosmological bound on neutrino masses is comparable 
and even stronger than direct terrestrial limits, presently at the eV order. A direct 
experimental limit on the neutrino mass comes from the experiments searching for 
kinematical effects on the spectrum of tritium (3 decay ^4) {m v < 2.2 eV at 95% 
CL), a bound expected to be improved by the KATRIN experiment to reach 0.35 
eV. A complementary bound can be put on the neutrino masses if they are Majorana 
particles, from the non-observation c of neutrinoless double beta decay 

The cosmological effect of neutrino masses is important because presently we have 
experimental evidences of flavour neutrino oscillat ions, w hich come from solar, atmo- 
spheric and reactor neutrino experiments (see e.g. BZEB t). They are sensitive to mass 
differences between the three neutrino mass eigenstates (mi,m 2 ,m 3 ). Approximate 
3cr ranges for the squared mass differences are Am 2 tm = 1.2 — 4.8 x 1CT 3 eV 2 and 
Am 2 un = 5.1 — 19 x 10~ 5 eV 2 , respectively^. 

These values are perfectly compatible with a hierarchical scenario where the light- 
est neutrino state has mi ~ 0, m 2 ~ (Am 2 un ) 1//2 and m 3 ~ (Am 2 ^) 1 / 2 (or with an 
inverted hierarchy where m 3 ~ m 2 ~ (Am^) 1 ' 2 , separated by the small Am 2 un ). 
The largest neutrino mass would then be of the order (Am 2 tm ) 1 / 2 ~ 0.05 eV. However, 
it is still possible that the three neutrino states are degenerate, with masses much 
larger than the differences. In such a case, J2i m i = m i + m i + m 3 — 3m . Cosmol- 
ogy is the most powerful tool for constraining the mass scale mo in the degenerate 
scenario. Applying the result from Eq. Q, we get m < 0.33 eV. 

At the current experimental level it is yet possible to reach the needed precision 
to test the non-degenerate neutrino schemes described above. However, future CMB 
observations by WMAP and PLANCK, among others, combined with LSS data from 
larger galaxy surveys will enhance the cosmological sensitivity to neutrino mass. The 

c It has been claimed a positiv e ev idence of neutrinoless double beta decay, criticized by other 
authors. For a discussion, see e.g. ESI. 



pioneering calculation in ref. found that the combination of future CMB and LSS 
data, and in particular the Sloan Digital Sky Survey (SDSS) will push the bound on 
vn v to approximately 0.3 eV. 

More recently, Hannestad has updated the forecast analysis and concluded 
that the PLANCK CMB experiment combined with the SDSS data could measure a 
neutrino mass of 0.12 eV at 95% CL, almost reaching the values in the hierarchical 
scenarios. With an order of magnitude larger survey volume than SDSS, the limit 
could be pushed down to 0.03 — 0.05 eV. 

There exist other cosmological probes of neutrino masses that could, in the next 
future, reach similar sensitivities. One can use the weak gravitational lensing of 
background galaxies by interveni ng m atter to probe the mass distribution of the 
universe. It has been estimated | 51 | 512 t a q -y e y sens itivity on m v could be 

reached^ depending on the equation of state of dark energy. Alternatively, the authors 
of ref. describe how to use the distortions of CMB temperature and polarization 
spectra caused by gravitation lensing to get a limit on m u . Their calculations show 
that again a ~ 0.1 eV value could be reach with PLANCK data. 

The cosmological limit on neutrino masses also applies for neutrino mass schemes 
that explain the results of LSND, a short-baseline accelerator experiment which gives 
an excess of v e events. These results constitute an additional evidence of flavour 
neutrino conversions than can be explained if Am^ND = 0-2 — 7 eV 2 . The minimal 
way to accommodate for LSND is to include a fourth sterile neutrino, and it was shown 
| 45 | 54 | only the so-called 3+1 scheme, where one heavy neutrino state of mass 
777,4 — (AmLgNj-,) 1 / 2 — J2 m u exists, in addition to the three flavour neutrinos, is still 
marginally allowed. The LSND results would be checked by the ongoing MiniBoone 
experiment . 

We have discussed before that the four neutrino scheme is disfavoured by BBN. 
In addition, shortly after the WMAP results several papers (see e.g. 
emphasized that the bound J2 m u < 0.7 eV was in conflict with the allowed windows 
for ?7i4 in a 3+1 scheme. We have discussed that this limit is probably too restrictive. 
It turns out that, in any case, this value or that in Eq. © can not be applied in 
the LSND case because it was calculated assuming three neutrinos. A calculation of 
the upper bound on m u including four neutrino states was done in ref. . It was 
found that the effect of an extra neutrino can partially compensate the effect of the 
neutrino mass, leading to a bound J2 m u < 1-38 eV (95% CL, N eS = 4). This value 
does not exclude yet the LSND windows with the lower m 4 values. 



5. Conclusions 



The existence of a relic neutrino sea in the universe allows us to study bounds 
on neutrino properties that are complementary to terrestrial experiments. We have 
discussed two different bounds that have been improved after the recent WMAP data, 



namely those on the effective number of relativistic neutrino species and on the sum 
of neutrino masses. Both are nice examples of the connection between cosmology 
and particle physics. The bound on AN e s restricts the existence of new relativistic 
particles or new decay processes beyond the Standard Model. The limit on J2 m u 
fixes the absolute neutrino mass scale, to which the current experimental evidences 
on flavour neutrino oscillations are not sensitive. 
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